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We present our recent studies of the pseudo-critical temperature, T c , of QCD using domain wall 
fermions. Domain wall fermions have the advantage that they preserve exact ST/ (2) chiral sym- 
metry at finite lattice spacing in the limit that L s — > oo. The RBC Collaboration has performed 
a set of dynamical calculations at L s = 32 and N t = 8 using the Iwasaki gauge action with two 
light quarks (m/a — 0.003) and one strange quark (m s a — 0.037). A clear signal for the crossover 
transition can be seen in the light chiral susceptibility, as well as in the Wilson line. However, 
at L s = 32, the residual chiral symmetry breaking is not yet fully under control. We also present 
preliminary results from the HotQCD Collaboration with N t = 8 and L s = 96, where the effects 
of the residual chiral symmetry breaking are reduced compared to L s = 32. 
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1. Introduction 

The location and nature of the QCD phase transition has been extensively studied using lattice 
techniques with various different fermion actions ^ ||, |3[ |J]. Recently, the most detailed studies 
of the transition temperature have been performed with different variants of the staggered fermion 
action @, B], which do not preserve the full chiral symmetry of QCD at finite lattice spacing. 
Domain Wall Fermions (DWF) on the other hand, allow the realization of exact chiral symmetry 
on the lattice, at the cost of introducing of an auxiliary fifth dimension. 

A study of the transition temperature with domain wall fermions was done for N t = 4,6 
However, it was found that at such coarse lattice spacings, the DWF formulation begins to break 
down, with unphysical effects that prevent the extraction of a reliable estimate of T c . 

In this work, we present a study by the RBC Collaboration of the pseudo-critical temperature, 
T c , with domain wall fermions with N t = 8 and fifth-dimensional extent L s = 32. It is hoped that, 
at the finer lattice spacings needed for N t = S, the large lattice artifacts that appear at N t =4,6 
are under better control. We also present preliminary results from the HotQCD Collaboration with 
domain wall fermions at N t = 8 with L s = 96. 

2. Simulation Details at L s = 32 

For our study we utilize the standard DWF action with an Iwasaki gauge action. The behavior 
of DWF at zero temperature has been extensively studied for this combination of actions; see refs. 
[Q, ^] for more details. 

The Rational Hybrid Monte Carlo (RHMC) algorithm [0, §], an exact algorithm that satisfies 
detailed balance and reversibility, is used to generate the gauge configurations. A three-level nested 
Omelyan integrator is used in the molecular dynamics evolution, with A = 0.22. The length of the 
molecular dynamics trajectories between Metropolis steps is chosen to be T = 1. The step size is 
tuned to achieve an acceptance rate of approximately 75%. A spatial volume of 16 is used for the 
finite temperature ensembles with N t = 8. For each value of the gauge coupling, we use a fixed 
value for the bare light and strange quark masses (mja = 0.003 and m s = 0.037). 

1200 molecular dynamics trajectories were also generated at j8 = 2.025 with a volume of 
16 3 x 32 and L s = 32, at the same quark masses used for the finite temperature ensembles. These 
configurations were used to calculate the static quark potential, as well as the meson spectrum. 

From the static quark potential, we obtain a value for the Sommer parameter, ro/a = 3.08(9). 
Using ro = 0.469(7) fm., this indicates a lattice scale of a~ ! rj 1.3 GeV at j8 = 2.025. The meson 
spectrum measurements, give a pion mass m % w 310 MeV, while the kaon mass is within 10% of 
the physical value. 

3. Finite Temperature Observables 

The observables that we use to probe the chiral properties at a given temperature are the light 
and strange quark chiral condensates ((i/h///) , (\jf\lf s )), and the disconnected part of the chiral sus- 
ceptibility (%i,Xs)- They are defined as: 

, > dlnZ 1 , . , s , 

^ = ^=m {Tr{<)) {3A) 

X q = ((W 9 ) 2 }-{W„} 2 (3-2) 
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Figure 1: On the left, (yv) f° r L s = 32.64.96. On the right, the disconnected chiral susceptibility for 
L s = 32,64,96. 

On all of our finite temperature configurations, we measure both the light and strange chiral con- 
densates every fifth trajectory, using 5 random sources per configuration to estimate Tr(M~ x ). 

Figure |l| shows the chiral condensate and the disconnected part of the chiral susceptibility, 
respectively. Examining the light and strange chiral condensate, it is difficult to precisely locate an 
inflection point, which is the signal for a thermal crossover. However, we can use the disconnected 
chiral susceptibility, a measure of the fluctuations in the chiral condensate. As seen in figure [j], 
there is a clear peak in the light disconnected susceptibility. The results for the chiral condensates 
and the associated susceptibilities are summarized in table [jj Fits to the peak region indicate that 
P c = 2.041(5). 

We also measure the observables that probe confinement, i.e. the Wilson line and its associated 
susceptibility. These are defined as: 

W^I^llj 1 ^ Xw = (W 2 ) - (W) 2 . (3.3) 

The results for the Wilson line and Wilson line susceptibility are also given in table [TJ. 
4. Residual Mass 

One of the primary drawbacks of the current calculation is the rather large residual chiral 
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Table 1: Summary of finite temperature observables, as well as the number of trajectories generated. 
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Figure 2: On the left, the dependence of m res on j3. On the right, the dependence of (y/y/) on L s at fixed 
input quark mass m q a = 0.003 at j3 = 2.0375. 

symmetry breaking for the parameters that we employ. This manifests itself in a value for the 
residual mass, m res which is larger than the input light quark mass, m/a = 0.003 over almost the 
entire range of parameters in our calculation. 

We have measured m res at j3 = 2.025 on the 16 3 x 32 zero temperature ensemble, which gives 
m res = .00665(8). We have also measured the residual mass on the finite temperature lattices. 
These measurements agree well with measurements on zero temperature lattices at nearby j8. 

Figure ^| shows how m res varies with j8. As we can see, m res has a strong, exponential depen- 
dence on p. While we have chosen the input light quark mass mi = 0.003 to be fixed at the different 
/3, the exponential dependence of m res means that the effective light quark mass, m q = mi + m res 
changes significantly in the crossover region, from m q « .0075 at j8 = 2.05 up to m q m .013 at 
P = 2.00. 

5. Chiral observables at varying L s 

The shifting of the quark mass with /3 results in a distortion of the shape of the susceptibility 
curves that we use to locate the crossover transition. In order to understand how this varying mass 
affects our results, we have measured the partially quenched chiral condensate with different L s 
and mi at various j3. In particular, we choose L s = 64 with the same input light quark, m/ = 
0.003, while for L s = 96, we vary the input quark mass so that the total effective quark mass m q 
approximately matches that for L s = 32 at the chosen value of j3. For one value of the gauge 
coupling (j6 = 2.0375), we measure y7y/ with many choices of valence L s and fixed input quark 
masses (mi,m s ) = (0.003,0.037). Table ||give the results of these measurements. Figure [l] shows 
the results with valence L s = 64 and L s = 96 in context with the L s = 32 results. 

From Figure [I], we see that holding the input quark mass fixed at mi = 0.003, while reducing 
m res by setting L s = 64 does not have a significant effect on the chiral condensate. On the other 
hand, with a larger input quark mass and L s = 96, the chiral condensate shifts appreciably. Figure 
^| shows the dependence of (xjfy) on L s at j8 = 2.0375. For small values of L s , there is a strong 
dependence, but the chiral condensate quickly plateaus to an approximately constant value for 
L s = 32,64, even though m res is still changing significantly in this region. 

In contrast to the chiral condensate, the disconnected part of the chiral susceptibility depends 
on the total effective quark mass, m q = mi +m res . As seen in figure [I], when the total quark mass 
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Table 2: Partially quenched measurements of (y/yf) at different L s , ft m/. 



is changed at L s = 64, the chiral susceptibility differs significantly from the measurements at L s = 
32. However, when we keep the total quark mass mi + m res fixed at L s = 96, the resulting chiral 
susceptibility agrees with L s = 32. Thus, while the chiral condensate is sensitive to the relative 
contributions from the input quark mass and the residual mass, the chiral susceptibility is a function 
of only the total quark mass, m q = mi + m res . 

6. Determining T c at L s = 32 

While the peak in the light chiral susceptibility is well-determined to be ft. = 2.041(5), there 
are several issues that need to be adressed in extracting a physical value for the pseudo-critical 
temperature, T c . 

Since m res changes so drastically as a function of j8, the chiral susceptibility curve is distorted 
by the changing light quark mass. Taking a simple ansatz %i ~ l/m q , we can adjust our data 
so that the bare quark mass is fixed as /3 varies. This adjustment shifts ft to stronger coupling, 
ft = 2.031(5). 

We have determined the lattice scale at j8 = 2.025, which differs slightly from ft. = 2.031(5). 
Using a simple interpolation between our result at j8 = 2.025 and the results at weaker coupling^], 
we obtain ro/a = 3.12(13) at j8 = j8 c . The results in ref. |Q] also indicate that chiral extrapolation 
and finite- volume effects add 4% to this value, giving ro/a = 3.25(18), where the error bar has 
been artificially inflated to include this 4% in the uncertainty. 

Since this calculation is done only at N t = 8 and at one set of quark masses, we cannot perform 
either the chiral or continuum extrapolation needed to obtain a value for T c at physical quark masses 
in the continuum. From ref. we can estimate the effect of the chiral extrapolation to the physical 
quark masses to be about 5%. Ref. §\§ also found that the effect of the continuum extrapolation is 
approximately 5%, although it is obtained using the p4 action. For our purposes, we estimate the 
error from the lack of continuum and chiral extrapolations to be 10%. 

Taking these errors into account, we obtain a value of T c tq = .406(23) (41), or T c = 17 1(10) (17) 
MeV, where the first error bar takes into account all the systematic errors outlined above except for 
the chiral and continuum extrapolation, which are reflected in the second error. 
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Table 3: Finite temperature observables for L s = 96 



7. Results at L s = 96 



The primary drawback of the RBC Collaboration's calculation just described is the rapidly 
changing residual mass as a function of /3 in the transition region. As we have discussed, this means 
that the effective quark masses in the strong coupling side of jS c are significantly larger than those 
on the weak coupling side. This has the effect of distorting the shape of the chiral susceptibility 
peak in the L s = 32 calculation, making the peak appear sharper and at weaker coupling than if we 
worked at fixed quark mass. 

The HotQCD Collaboration's calculation seeks to address this flaw by working at L s = 96. 
Utilizing L s = 96 reduces m res approximately by a factor of 3. It also allows us to choose the input 
quark masses (mia and m s a) so that the total effective quark mass ((m; + m res )a and (m s + m res )d) 
are the same at each value of j8 that is used. For this calculation, the light quark mass is chosen 
to be 1/6 the strange quark mass at each value of j8, where the strange quark is chosen to be 
approximately physical. The corresponding bare quark masses are given in table ||. Otherwise, the 
lattice actions, the spatial volume, and the molecular dynamics algorithm used are all identical to 
those used at L s = 32. 

Table |5| gives preliminary results for the chiral condensates, wilson line, and their susceptibil- 
ities. These results are also presented in figure |3[ As expected, working at fixed quark mass results 
in a chiral susceptibility where the peak is broader and much harder to resolve. In addition, there 
are some indications that /3 C is at stronger coupling at L s = 96 compared to L s = 32, as expected, 
but the statistical error makes this far from a certain conclusion. 



HotQCD 
Preliminary 
(L=96) 



g-o Ls = 96 pbpsus 
B-a Ls = 32 pbp_sus 




o-e> Ls = 96 pbp - 
h-oLs = 32 pbp 
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Figure 3: On the left, a comparison between L s = 32 and L s = 96 of the light quark chiral condensate and 
chiral susceptibility. On the right, the same comparison for the Wilson line and its susceptibility. 
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8. Conclusions and Outlook 

We have presented two studies of the critical region of finite temperature QCD using domain 
wall fermions at N t = 8. The calculation at L s = 32 by the RBC uses observables related to chiral 
symmetry (i.e. chiral condensate and disconnected chiral susceptibility), to calculate the crossover 
transition temperature, giving the result T c r = .406(23)(41). or T c = 171(10)(17) MeV. 

The second calculation by the HotQCD calculation improves upon the RBC calculation by 
using L s = 96 to further reduce the residual chiral symmetry breaking, while tuning the input quark 
masses so that the total bare quark mass (including m res ) is fixed at each different value of /3. 
However, preliminary results show that the peak is less sharply resolved compared to L s = 32, so 
no reliable estimate for /3 C can yet be obtained. 

Data at a few additional j3 at both stronger and weaker coupling are needed to better resolve 
the shoulders of the chiral susceptibility peak (if indeed such a peak exists). In addition, a zero 
temperature calculation to determine the lattice scale needs to be done in order to obtain a physical 
value for T c in MeV. 
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